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Аннотация 

Введение. Пульсации рабочей среды в гидросистемах станков, системах топливоподачи авиационных 
двигателей и подачи жидкостного ракетного двигателя являются одним из главных факторов возникновения 
вибрации агрегатов, приводящей к усталостному разрушению корпусов элементов, трубопроводов и выходу из 
строя элементов насосов. Целью работы являлись реализация метода определения динамических характеристик 
объёмного насоса с использованием специальных стендовых систем и сравнение результатов расчёта с 
экспериментальными данными. Описаны этапы расчёта динамических характеристик объёмного насоса с 
учётом разработанных ранее специальных стендовых систем на примере шестерённого насоса внешнего 
зацепления производительностью 14 см/об. Показана реализация метода вариации нагрузок В. П. Шорина с 
использованием разработанных авторами специальных стендовых систем с заранее определёнными 
динамическими характеристиками. Излагаются основные этапы методики определения динамических 
характеристик шестеренного насоса. 

Материалы и методы. В работе применены методы спектрального анализа пульсирующего давления. 
Пульсации расхода жидкости на выходе из насоса определялись с помощью импедансного метода, метода 
вариации нагрузок и специальных стендовых систем. 

Результаты исследования. В работе реализована методика определения динамических характеристик 
шестеренного насоса в диапазоне оборотов приводного вала 500-2500 об/мин для четырёх гармонических 
составляющих спектра колебаний в широком диапазоне динамических нагрузок (от инерционного до 
емкостного характера). Представлен анализ стендовых систем, приводящих к расчёту динамических 
характеристик насоса с минимальной погрешностью на основании условия согласования динамической 
нагрузки и источника колебаний. Проведена проверка разработанного подхода к оценке динамических 
характеристик насоса посредством сравнения расчётных и экспериментальных данных пульсаций давления в 
стендовых системах с дросселем, полостью и протяжённым трубопроводом на выходе насоса. 

Обсуждение и заключения. Реализован метод определения динамических характеристик объёмного насоса с 
использованием разработанных авторами специальных стендовых систем. Результаты исследований 
показывают, что исследуемый шестерённый насос может рассматриваться как независимый источник 
колебаний расхода, для которого отклонение собственных динамических характеристик от средних значений не 
превышает 10 % по первой гармонической составляющей. 


Ключевые слова: объёмный насос, метод вариации нагрузок, стендовые системы, источник колебаний расхода, 
динамические характеристики насоса, независимый источник колебаний, пульсации расхода, импеданс. 


Финансирование: работа второго автора выполнена по теме № Е555-2020-0015 в рамках выполнения 
госзадания Минобрнауки России. 


© Санчугов В. И., Рекадзе П. Д., 2022 


Санчугов В. И. и др. Определение динамических характеристик шестерённого насоса по методу вариации нагрузок 


Для цитирования: Санчугов, В. И. Определение динамических характеристик шестерённого насоса по методу 
вариации нагрузок с использованием специальных стендовых систем / В. И. Санчугов, П. Д. Рекадзе // Ад4уапсеа 
Епотеегте Везеагсв. — 2022. — Т. 22, №2. — С. 130-141. №&рз://А01.0г®/10.23947/2687-1653-2022-22-2-130-141 


Опвта[ агисе 
Реегитайоп о Фе дупапис свагасег1$9с$ оРа геаг рипр Бу 1е 1оа4 
уачаНоп ше@од ито 5реса! Бепсв $у$ет$ 
Узену Г. Запсви?оу © Рауе! О. Кекад7е Ох 


Затага МаНопа] Кезеагсн ОшуегзИу, 34, МозКоузКоуе 5Воззе,батага, Кизз1ап Недегайоп 


Ч ВеКаа7е1993 @ ятай.сот 


АЪ5етгасе 

Гитодисйоп. Зоте оЁ Фе та]ог Расюг$ ш Фе оссштепсе о{ угайоп оЁ ипи$, сало Гайеие аПиге оЁ Фе Вои$115$ 
ог @етеп($, ргреПпез, ап РаИите оЁ ритр @етеп, аге р 5айопз оЁ Фе мойте теапит ш е вудгаиНс зузета$ оЁ 
тасбше 10015, Рае! Еее зууеплз оЁ ашсгай епошез ап4 Паи1А-ргореПапе епоше зарр1у. ТЫ5$ зеа4у аппе4 аЁ ше 
пиретешаНоп оЁ а шефо4 Гог аеегиште Фе Чупапис спагасетзс$ оЁР а уотейлс ритр изше зресла! Бепсв 
зузет$, ап4 Ше сотра5оп оЁ Ше са1сшШайоп гези5 ю Ше ехрегипета| Даа. ТБе $%асез оЁ сайсшШайпя фе дупапис 
спагасет1$Яс$ оЁ а уоштеелс ритр \ете дезспЪеа, {акте шю ассоцпЕ Ше рге-4еуе]оре4 зрес1а1 БепсН зубет$ оп 
Ще ехатре оЁ ап ех(егпа! эеаг ритр \МИВ а сарасНу оЁ 14 сш”/еу. Те паретепканоп оЁ У. Р. $Вони'з 1оаа 
уапаноп тешо4 изшз зреста! БепсВ 5узепл$ 4еуеоре4 Бу Фе аиог$ \ИВ ргедеегилте4 дупапис сБагасетзНс$ \аз 
зВо\п. ТВе тат ${асез оЁ ше теодо]оэу Гог ае®егитите Фе дупапис сВагасетзИс$ оЁ а сеаг рштр \ете АезспЪеч. 
Миепта[5 апа Мейо4$. Мефод$ о{ зресга| апа1уз1$ о ршзаНпе ргеззиге \еге аррПеЯ ш Фе \отк. Ри15айопз 
ог Иша По\ а Фе ритр оц е! уеге ащетттеЯ ито Ше ппредапсе те!фо4, ше те'о4 оЁ 1оа4 уапайоп, ап4 
зресла1 Бепсв зузетт$. 

Вези5. Тре рарег ппр!етепие4 а {есбшаче ог 4еегтиипе Фе дупапис свагасет5с$ оРа зеаг ритр ш Фе 4пуе 
зваЁ зрее гапзе оё 500-2500 грит Гог юг Вагтошс сотропепз оЁ Фе угайоп зресигат ш а \14е гапсе оЁ аупапис 
1юа4$ (Нот шегНа| ю сарасуе пашге). ТВе БепсН зузетз уле!Ате Ше са1саНоп оЁ фе дупапис сВагасет$с$ оЁ 
фе ритр \1 а шшипит еггог Базе оп ше соп@! оп оЁ тасЫше Фе дупапис 1оа4 ап Ше зоигсе оЁ угаНопз \еге 
апа]у7е4. ТЬе 4еуе!оре4 арргоасв 0 Ше еуашаНоп оЁ Ше Аупапис сВагас{е1$Ис$ оЁ Ше ритр \аз уемНе4 Шгоч>В 
сотрание Фе са]сае4 ап ехреттепиа| даёа оЁ ргеззиге ршайопз ш @е БепсН зузет$ м своКе, сауПу ап4 ап 
ежепае4 р:ре!пе а{ Фе рштр оиЙей. 

Оёсизяют апа Сопсш51оп5. ТВе тефо4 Гог деегииите Фе 4упапис сБагасет$Яс$ оЁ а уоатешлс ритр \аз 
пир!етеше4 изшс зреса| БепсВ зузеил5 4еуе]оре4 Бу Ше аиогз. ТВе гезеагсВ гези{$ зВо\/ фа{ Ше зеаг рштр ипаег 
з@а4у сап Бе сопз1еге4 аз ап ш4ереп4епе зоигсе оЁ Но\ НисваНоп$, ог \ЮсЬ Ше Аеуайоп оЁ 15 омп упапис 
свагас(ег1$с$ Гот Ше ауегасе уащез 4оез по ехсееа 10% Гог Фе НгзЕ Вагтошс сотропепе. 


Кеумога$: уоитевлс ритр, юа@ уапаноп тефо4, БепсВ зузетз, зоигсе оЁ Йом НасшаНоп$, ритр Чупапис 
спагасет1зИс$, ш4ереп4еп{ зопгсе оЁ НисаНопз, Но\ ри5аНопз, ппредапсе. 


Еипаше и!огтаНоп: Ше гезеагсВ оЁ Фе зесоп4 айог \'а$ сагле опё оп Фете по. Е555-2020-0015 уиш е 
Нате оЁ Фе гоуегитеп! (азК оЁ Ше Миизку оЁ Едисайоп апа Заепсе оЁ Фе Киз1ап Не4дегаНоп. 

Ког сНабоп: У. Г. Запсвигоу, Р. О. ВеКадте. Реегиитаноп о{ Фе Аупапис сБагасетзЯсз оЁа сеаг ритр Бу Ше 1оа4 
уапаноп тешфо4 изше зрес1а! БепсЬ зузетз. А4уапседа Епошеетто КезеагсН, 2022, уо1. 22, по. 2, рр. 130—141. 
№ рз://401.0г2/10.23947/2687-1653-2022-22-2-130-141 


Введение. Основные недостатки насосов заключаются в сопровождаемом их работу акустическом шуме, 
вибрации корпуса, присоединенных трубопроводов и возникающих вследствие этого разрушениях конструкции. 

Традиционный путь совершенствования конструкций насосов, используемый в машиностроении и 
заключающийся в сглаживании неравномерности течения жидкости, затруднен отсутствием средств прямого 
измерения расхода жидкости на выходе насоса. Это ограничивает любую экспериментальную оценку 
эффективности конструктивно-доводочных работ по насосам. Кроме того, это ограничивает аналитическое 
описание насоса как источника колебаний, необходимое разработчикам гидросистем для обеспечения 
требуемого динамического качества самих систем. 

По этой причине исследования, связанные с моделированием динамических процессов в гидромеханических 
системах [1,2] и применением новых подходов к решению задач динамики гидравлических систем [3, 4] на 
основе физических принципов, имеют актуальность. 
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Материалы и методы. Представлена реализация метода вариации нагрузок В.П. Шорина [4], 
заключающаяся в том, что определяются пульсации давления за насосом, взаимодействующим с несколькими 
стендовыми системами, динамические характеристики которых рассчитаны заранее. Далее обработанные 
амплитудно-частотные характеристики давления с использованием импеданса стендовых систем 
пересчитываются в собственные динамические характеристики насоса — пульсационную производительность 
и импеданс. Развитие подхода состоит в том, что предложены специальные стендовые системы, учитывающие 
соединительные штуцеры, переходники и внутренние каналы насоса, что позволяет точнее определять 
динамические характеристики насосов, а также осуществлять проверку применимости используемой модели 
насоса в виде эквивалентного источника колебаний. 

Результаты исследования 

Методика определения динамических характеристик шестеренного насоса 

В Самарском университете разработана расчетно-экспериментальная методика определения динамических 
характеристик шестеренного насоса (пульсационной производительности и импеданса), которая базируется на 
исследованиях Л. Тевенена, Э. Нортона и В. П. Шорина [5, 6]. Разработанная методика включает в себя 
следующие основные этапы: создание моделей и расчет динамических характеристик специальных стендовых 
систем; экспериментальное определение пульсаций давления на выходе насоса; расчет спектра возбуждаемых 
колебаний в различных стендовых системах; определение аппроксимирующих зависимостей для отдельных 
гармонических составляющих; расчёт собственных динамических характеристик насоса (пульсационной 
производительности и импеданса); оценка независимости (стабильности) динамических характеристик насоса 
от характеристик стендовых систем; проверка применимости используемой модели насоса в виде 
эквивалентного источника колебаний расхода (или давления). Далее приведен пример реализации 
предложенной методики с использованием шестерённого насоса внешнего зацепления удельной 
производительностью 14 см“/об. 

Модели и расчет динамических характеристик специальных стендовых систем 

Рассмотрены базовые варианты стендовых систем на выходе жидкости при испытаниях насосов. В 
сосредоточенных параметрах стендовая система может реализовывать активную расходно-зависимую, 
инерционную и упругую нагрузки. В распределённых параметрах стендовая система может быть реализована в 
виде протяжённого цилиндрического трубопровода, реализующего активную частотно-независимую 
неотражающую нагрузку. 

Модели насоса в виде эквивалентных источников колебаний расхода (ЭИКР) или давления (ЭИКД) 
определяются с использованием источника постоянного расхода О,, идеального источника колебаний давления 


р (или расхода 4) и внутреннего импеданса 7, [4]. На рис.1 представлены модели насоса в виде 


эквивалентного источника колебаний давления (и расхода) с присоединенной системой. 


стендовоя 
НаСОС | система насос | система 
а) 6) 


Рис. 1. Модели насоса в виде эквивалентного источника колебаний: а — давления; 6— расхода с присоединенной 


системой: О, 5 Б, = средний расход и давление; р„,4и— собственная пульсационная производительность насоса по 


давлению и расходу, 7, Иа — импеданс насоса и стендовой системы (адаптировано из [4]) 


и? 


При этом задачей выбора стендовой системы является обеспечение таких параметров системы, при которых 
реализуются наилучшие условия высокой точности измерений. Задачу выбора параметров присоединенной 
системы необходимо проводить при изменении импеданса в диапазоне от -© до +00. 

Принципиальные схемы стендовых систем представлены на рис. 2. 
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Рис. 2. Принципиальные схемы стендовых систем для определения собственных динамических характеристик насосов: 
а — сактивной дроссельной нагрузкой; 6 — с активной неотражающей нагрузкой; в — емкостной нагрузкой: 
1 — испытуемый насос; 2 — нагрузочный дроссель; 3 — расходомер; 4 — датчик давления; 5, 6 — соединительные 
штуцеры; 7 — магистраль; 8 — протяжённый трубопровод; 9 — полость 


Расчётные модели импедансов специальных стендовых систем: 
— с активной дроссельной нагрузкой (дроссель): 
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— с активной распределенной нагрузкой (протяжённый трубопровод): 
И & ты _— 2 _ 


тр 


— с емкостной нагрузкой (полость): 


2 (3) 
те) 1-1 
С тех (О, Л ) = (ке? (0. ) -+ Кей ив (0, ) + Кей. (0. ) + р (0. )} + ЕЕ ы 
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В формулах обозначены: 7.» (0,) —Щ импеданс дросселя; Кей’ (©) —щ действительная часть импеданса 


трубопровода, включающего соединительную магистраль, штуцеры, внутренние каналы в насосе и дросселе; 
Кей ‚Кей _ ‚ Кей —Й действительная часть импеданса штуцеров и внутреннего канала в 
`лам о `турб о 'м.с. о 


насосе и полости при ламинарном, турбулентном течении жидкости за насосом и местных сопротивлениях 
соответственно; {, — частота К-й гармонической составляющей спектра колебаний за насосом; р — 


Плотность рабочей среды; 1; —Щ длина трубопровода, включающая длину магистрали 6, соединительных 
штуцеров 5и7и внутренних каналов насоса и дросселя (не показаны на рис. 7); ой — площадь поперечного 


сечения трубопровода; Г, — «инерционность» соединительных штуцеров и внутреннего канала в насосе и 
полости; С — упругость полости; К — номер гармонической составляющей спектра колебаний; | — мнимая 


единица (] = -1 ). 
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С использованием стендовой системы с дросселем на выходе насоса реализованы три инерционно-активные 
динамические нагрузки, в которых изменялось активное сопротивление изменением площади перекрытия 
дросселя Но: 

— для стендовой системы №1 площадь перекрытия Ео!=9,290х 10° м2, 

— для стендовой системы №2 площадь перекрытия Ео›=7,299х 10° м7, 

— для стендовой системы №3 площадь перекрытия Е‹з=2,247х 103 м?. 

С использованием стендовой системы с полостью на выходе насоса объёмом 2 л реализована инерционно- 
упругая динамическая нагрузка. 

С использованием стендовой системы с протяжённым трубопроводом длиной 106 м реализована активная 
частотно-независимая нагрузка. 

По моделям стендовых систем по формулам (1)—(3) рассчитаны зависимости импеданса (модуля и фазы) от 
частоты. В таблице 1 для 1-й гармоники приведены функции регрессии модуля импеданса стендовой системы 
7(Р и фаз В;(Е) для анализируемых стендовых систем. 


Таблица 1 
Функции регрессии динамических характеристик стендовых систем по 1-й гармонике 
№ Модуль импеданса 
стендовой Тип стендовой стендовой системы Фаза импеданса стендовой системы 
системы по системы кг В; (В, рад 
Я: 
рисунку 2 м‘с 
а— №1 1,412-107+ 2,201-107 7,066-10# + 0,199 
С дросселем на - 9 т 
а— №2 НЕОО-НАСОСЬ 1,610-10 + 4,258:10 6,906:10"Е+0,09039 
а— №3 4,980-107Ё + 1,503-10? 2,723-103#+ 0,251 
С полостью на 3,308-10“Ё — 1,416.10 + 4,234-10° Г +1,495.10? 
б- №4 9 атс ( - с) 
выходе насоса 1,873.10 1,430.10} +2,930.10 
С протяжённым 
в- №5 трубопроводом на 1,897-1010 0 
выходе насоса 


Экспериментальное определение пульсации давления на выходе насоса 

Исходной информацией для определения динамических характеристик насоса являлось измерение давления 
на выходе насоса. Пример осциллограммы и спектра пульсаций давления на выходе из насоса в стендовой 
системе с активной распределенной нагрузкой представлены на рис. 3. Режим работы насоса при этом: среднее 
давление РД, = 23,8 МПа ‚ частота вращения приводного вала п =100006 / мин. 


24,30 уе =-—--------- ртахЕ24,26 МПа-------------- и и 
алз 1 АНН Мм 
23,95 | 
> 
9“ 23,8 
Е 
ы | | | 
2 23,43 обои о бои К РО С о | пои у и О о У Я © о 
23,25 = 1 ооо | ме ры Зоо И | у ТИ 
23,08 о ЗВАНИИ: БЫ ВЕ ЗЕ Бип: 2305 МП | 
Т ы Т т т ы Т т с ы 
0,03 0,04 0,05 0,06 0,07 0,08 
Время поворот ротора, с 
| т т Т Т Т | 
0 60 120 180 240 300 360 
Угол поворота ротора, град 
а) 
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Амплитуда пульсаций давления, МПа 


167 334 501 668 835 1002 1169 1336 1503 1670 1837 
Частота составляющих спектра колебаний, Гц 
6) 


Рис. 3. Давление на выходе насоса: а — осциллограмма; 6 — спектрограмма 


Амплитудно-частотная характеристика пульсаций давления на выходе насоса в диапазоне оборотов 
приводного вала насоса п=500-1250 об/мин представлена на рис. 4. Сплошной линией нанесены функции 
регрессии, построенные в среде Мсгобой Ехсе! при помощи инструмента «Линия тренда». 


Ь: 

3 0,6 

> ® | гармоника 
Е 0,5 ©2 гармоника 
: 0.4 © 3 гармоника 
я ® 1 гармоника 
Е 

= 0,3 

з 

: 

в. 02 

з 

5 

2 0,1 

Е 

Е 

з‹ `0 


83 166 249 332 415 498 581 664 747 830 


Частота составляющих спектра колебаний 13, Гц 


0 500 1000 1500 
Частота вращения приводного вала п, об/мин 


0 7,3 14,6 21,6 


Средний расход жидкости насоса О, л/мин 


Рис. 4. Амплитудно-частотная характеристика пульсаций давления на выходе насоса по отдельным гармоническим 
составляющим (точками обозначены данные на: ®@ — первой, ® — второй, ® — третьей, ® — четвёртой гармониках, 
линиями обозначены функции регрессии) 


Анализ осциллограммы, спектрограммы и амплитудно-частотной характеристики пульсаций давления 
показал следующее: тип процесса полигармонический установившийся; число регистрируемых гармоник — 8, 
анализируемых' — 4; диапазон амплитуд колебаний давления 0,28—0,52 МПа; частотный диапазон приводного 
вала 500-1250 об/мин, анализируемых гармонических составляющих спектра 83-833 Гц. 

Для гармонического анализа периодических сигналов использовалось разложение в ряд Фурье [7]. Для 
разложения в ряд Фурье использованы участки осциллограммы давления, измеряемые относительно опорного 
сигнала на роторной частоте, т.е. за полный оборот ротора извлекался один участок осциллограммы. Опорный 
сигнал регистрировался при фиксированном положении шестерни, выставляемом при начале записи 
осциллограммы давления для пяти стендовых систем. Таким образом определялись амплитудно-частотные 


(АЧХ) А, ( @) и фазо-частотные (ФЧХ) ф, (<) характеристики. 


1 В работе принято, что крайняя высшая анализируемая гармоника по амплитуде не должна быть более 5 % от амплитуды первой 
гармоники. 
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Расчет спектра возбуждаемых колебаний в различных стендовых системах и определение 
аппроксимирующих зависимостей для отдельных гармонических составляющих колебаний давления 

На основе экспериментального определения давления в пяти стендовых системах: с дросселем (№1-3), с 
полостью (№4) и протяжённым трубопроводом (№5) на выходе насоса рассчитаны спектры возбуждаемых 
колебаний. В таблице 2 приведены функции регрессии амплитуд А,(Р) и начальных фаз пульсаций давления для 
1-й гармоники. 


Таблица 2 
Функции регрессии амплитуд и начальных фаз пульсаций давления 
№ стендовой . Амплитуда пульсаций давления Начальная фаза р 
Тип стендовой ы Частотный 
системы по кг пульсаций давления 
системы Ар (в, —— диапазон 
рисунку 2 мс ф;: (В, рад 
а— №1 1,003-107% —0,543 
С дросселем на Е 
а— №2 выходе насоса 1,243.10 Е —0,280 
а— №3 6,129.10 —5,941 
С 
6 №4 о 1,343-10°°Р — 4,298-104#+ 3,913-107 -5,373 83—417 
выходе насоса 
С протяжённым 
в №5 трубопроводом на 2,694-107% —0,658 
выходе насоса 


Расчёт динамических характеристик насоса 


Для расчета динамических характеристик насоса (импеданса и пульсационной производительности 
источника) необходимо провести эксперимент как минимум с двумя динамическими нагрузками”. При этом 
импеданс насоса по уточнённой авторами модели эквивалентного источника колебаний расхода с 


использованием аналитических зависимостей р, (/ ) и 7, (1) из таблиц 1 и 2 рассчитывается по формуле: 


1 
=2, а. (4) 

и И 

ры 

р. 2, 

Г 2 2 

где 7, — импеданс насоса (7, =|7, е*) — модуль импеданс насоса (|2, =\/Кейи +Пий и ), Вей „1т7, — 
действительная и мнимая часть импеданс насоса, ф — аргумент импеданса насоса, 7,,/, — входной 
импеданс динамических нагрузок № 1-2 (= е^, 2, =|2.|е*), 8,8, — аргументы имнедансов 
динамических нагрузок № 1-2, р,р, —щ пульсации давления за насосом в стендовой системе № 1-2 ( 
р = Ане" ‚Р› =А зе ), Ат, Ар» — амплитуды пульсаций расхода за насосом в стендовой системе № 1-2, 


ф,ф, — начальные фазы пульсаций давления динамических нагрузок № 1-2. 

Затем с использованием аналитических зависимостей из таблиц 1,2 определяется пульсационная 
производительность насоса по уточнённой авторами модели эквивалентного источника колебаний (переменная 
составляющая объёмного расхода а, и давления р, ): 


Ее 1-2 
о (5) 
“ 2, | и И 
Р Р> 2 
где 4,‚р, — переменная составляющая объёмного расхода (а, = Ане) и давления (р, = Аше” ”_) насоса, 
А„.Ам — амплитуда пульсаций расхода и давления насоса, ф,ф, — аргумент переменной составляющей 


объёмного расхода и давления насоса. 
Из рассчитанных характеристик по модели ЭИКР (а, ( Г) и ЭИКД (р, ( У )) по формуле (5) только одна 


будет отличаться стабильностью, т.к. эквивалентный источник колебаний имеет независимую динамическую 
характеристику только по одному из параметров (пульсирующий расход или давление). 


2 И ы 
Для достаточного объёма выборки необходимо использование хотя бы четырёх динамических нагрузок 2,, формирующих шесть 


неповторяющихся (в том числе 1 = ] ) комбинаций: бо, бл ба» Из, бд, Иа. 
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Расчет импеданса 7, и переменной составляющей объёмного расхода 4, насоса по формулам (4)-—(5) 


проведён для пяти динамических нагрузок, т.е. с использованием десяти различных комбинаций этих нагрузок: 
«1_2», «1_З», «2_З», «14», «2_4», «З_4», «1_5», «2_5», «З_5», «4_5» (см. таблицу 1). 

Для дальнейшего расчёта и анализа из десяти кривых амплитуд пульсаций расхода получена характеристика 
амплитуд пульсаций расхода насоса А„. Значения амплитуд А„ (на каждой частоте) рассчитывались по 


формуле из источника [8]: 


А„=А„ АА, 
где А„ — средняя амплитуда пульсирующего расхода, ДА„, — доверительный интервал амплитуды пульсаций 
19 
расхода насоса (ДА =-——), о — среднеквадратическое отклонение амплитуды пульсирующего расхода от 


Ут 
среднего значения, и — число комбинаций динамических нагрузок (т = 10), { — коэффициент Стъюдента ( 


1=2,262 , рассчитан для доверительной вероятности 0,95). 

Доверительный интервал амплитуд ДА„ формируется из погрешностей (инструментальной и методической) 
и степени зависимости динамической характеристики насоса А, (} ) от стендовой системы. Инструментальная 
погрешность состоит из погрешности определения давления и фаз, методическая погрешность состоит из 
погрешности неучтённого влияния динамических свойств приводной механической системы. 

Фазы пульсаций расхода ф, и их отклонения рассчитывались аналогичным образом. При расчёте фаз 
комбинации «1_2» и «1! _5» исключены по причине вырождения переменной составляющей объёмного расхода 
насоса 4, при использовании близких значений характеристик на одной частоте: ф =фир =р,. 

Результаты расчёта собственной динамической характеристики насоса в виде пульсационной 
производительности и импеданса представлены на рис. 5. Фазы для удобства представлена в градусах. 
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Частота гармонических составляющих спектра колебаний, Гц 
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2 
= 156840 
м 
= 14Е+10 
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о 
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Частота гармонических составляющих спектра колебаний, Гц 
г) 
Рис. 5. Динамические характеристики насоса: а — амплитуда пульсаций расхода; 6 — фаза пульсаций расхода; 


в — модуль импеданса; г — фаза импеданса (сплошными линиями обозначены данные по: ча» — первой; 
чннинию — ВТОРОЙ; чине —— Третьей; чаю —— четвёртой гармониках) 


Амплитуда пульсаций расхода на 1-й гармонике имеет монотонно возрастающий характер. С 
увеличением номера гармоники наклон кривой амплитуд пульсаций расхода снижается до практически 
постоянного значения к 4-Й гармонике. Значения амплитуд А, на высших гармониках на порядок меньше, 


чем на 1-й гармонике. 
Модуль импеданса насоса на 1-Й гармонике имеет монотонно возрастающий характер. 
Для оценки реализованных динамических нагрузок использован критерий согласования модуля импеданса 


|2, источника колебаний и входного модуля импеданса динамической нагрузки |2 [99, 10]. Модули 


ст.с 


импедансов должны быть сравнимы (|7,|=| 


ст.с 


), чтобы количественные значения динамической 


характеристики насоса определялись с минимальной погрешностью. Для сравнения использовано отношение 


и и 


7 й При =1 динамическая нагрузка считается согласованной с источником. 


ст.с. ст.с. 


Ввиду различия форм кривых модулей импеданса насоса и динамических нагрузок в диапазоне частот 83— 


и 


417 Гц рационально использовать среднее отношение ‚› Которое позволит учитывать вклад большей 


Ст. С- 1 и Тнах 


части частотного диапазона. Т.е. отношение по 1-й гармонике составит: 
нах 

— для активной дроссельной нагрузки №1 — 1,16; 

— для активной дроссельной нагрузки №2 — 0,8; 

— для активной дроссельной нагрузки №3 — 2,65; 


— для емкостной нагрузки №4 — 4,64; 
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— для активной распределенной нагрузки №5 — 0,35. 

Т. е. динамические нагрузки стендовых систем №1,2 являются наиболее согласованными с источником 
колебаний. Таким образом, динамическая характеристика насоса, рассчитанная с использованием 
динамических нагрузок стендовых систем №1,2 («1 3», «23», «1 4», «24», «|5», «2 5») определена с 
минимальной погрешностью. 


Оценка стабильности динамических характеристик насоса от характеристик стендовых систем 


Стабильность динамических характеристик Ам и ф, ри, 


и ф,) насоса оценивалась с использованием 


инструментов математической статистики [8]. Так стабильность характеристики Ам оценена через расчёт 


отношения доверительного интервала [АА на частоте Тк среднему значению характеристики | на этой 
3; 
частоте Полученное относительное отклонение А, не должно превышать заданное отклонение Азот : 
ДА 
— Чи 
д. т А — Азот ? 
4иМе 
где Ашме среднее медианное значение амплитуды пульсаций расхода насоса; Аа задаваемое 


относительное отклонение амплитуды пульсаций расхода насоса от среднего значения. 


При относительном отклонении Д, —>0 (когда [АА 


—0) стабильность характеристики А„ максимальна. 


Исходя из опыта других исследователей, предпочтительно выбирать величину А, в пределах 5-30 % [11-17]. 


допа 


Для одночисловой оценки динамических характеристик насоса рационально использовать величину 
относительного отклонения. Расчёт такой величины возможен по оценке среднего медианного значения 


ДА, 
относительных отклонений в диапазоне частот каждой гармоники ( 1... Ли»): Ламе = = . Величину 
А. Лт --“Лтах 
АА х 
Ам рационально использовать, если — < 0,3 ‚ т.е. когда доверительный интервал ДА» не превышает 


44 У, 
30 % от средней амплитуды пульсирующего А„. 


Результаты оценки стабильности динамических характеристик насоса от характеристик стендовых систем“ 
составляют: по 1-й гармонике — 2%, по 2-й — 8%, по 3-й — 15%, по 4-й — 48%. Результаты расчёта 
динамических характеристик насоса показали, что по четырём гармоникам в диапазоне частот 500-2500 об/мин 
исследуемый шестерённый насос целесообразно рассматривать по модели эквивалентного источника 
колебаний расхода. 

Проверка применимости используемой модели насоса в виде эквивалентного источника колебаний расхода 

Для проверки разработанного подхода оценки динамических характеристик насоса было выполнено 
расчётное определение пульсаций давления в стендовых системах с дросселем, полостью и протяжённым 
трубопроводом на выходе насоса. С использованием найденных динамических характеристик насоса 
рассчитаны амплитуды пульсаций давления в анализируемых стендовых системах, которые сравнивались с 


экспериментальными. Ниже представлены результаты расчёта амплитуд А, ( 1) в стендовой системе с 
протяжённым трубопроводом. Амплитуды пульсаций давления р, рассчитаны в соответствии с моделью ЭИКР 
по формуле: 
4,:7,-7, (6) 
7,47, 
Результаты расчёта амплитуд пульсаций давления за насосом приведены на рис. 6. Пунктирной линией 
обозначены амплитуды, полученные экспериментально, а сплошной линией — рассчитанные по формуле (6). 


р; = 


ь Стабильность других собственных характеристик насоса оценивается аналогично. 
ь Необходимо учитывать методическую погрешность подхода и инструментальную погрешность определения давления. 
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Рис. 6. Амплитудно-частотная характеристика давления для стендовой системы с протяжённым трубопроводом длиной 
[=106 м (сплошными линиями обозначены расчётные данные по: ча» — первой; чение — второй; 
чниииинию — ТрстьсйЙ; чин — четвёртой гармониках; пунктирными линиями обозначены экспериментальные данные) 


Оценка сходимости результатов проверочного расчёта проведена по формуле отклонения разницы 
экспериментальных и расчётных значений от экспериментальных данных. Результаты проверочного расчёта для 
стендовых систем с дросселем, полостью и протяжённым трубопроводом на выходе насоса показали, что по 1-й 
гармонике относительное отклонение от экспериментальных данных для рассмотренных стендовых систем не 
превышает 10%, минимальное относительное отклонение характерно стендовой системе с протяжённым 
трубопроводом на выходе насоса — 2%. Амплитуды пульсаций давления 2-й-4-й гармоник находятся в пределах 
коридора, рассчитанного по формуле А, =А,-А, (где А, — средняя амплитуда пульсирующего давления). 

Обсуждение и заключения. Реализована методика расчёта динамических характеристик насоса в виде 
эквивалентного источника колебаний расхода жидкости на примере шестерённого насоса производительностью 
14 см*/об с использованием разработанных авторами специальных стендовых систем и уточнённой модели 
импеданса и пульсационной производительности по расходу. Показано, что исследуемый шестерённый насос 
внешнего зацепления рационально рассматривать по модели эквивалентного источника колебаний расхода. 
Проведена проверка предлагаемого подхода в виде расчёта амплитуд пульсаций давления в стендовых 
системах и сравнения с экспериментальными значениями. Отклонение расчётных амплитуд пульсаций 
давления от экспериментальных значений не превысило 10 %. 

Стоит отметить, что за рамками обзора осталось влияние на динамические характеристики шестерённого 
насоса приводной механической системы и износа деталей шестерённого насоса. На основе проведённой 
работы можно сформулировать направления развития предлагаемого подхода: 

— оценка динамических характеристик шестерённого насоса в течение всего срока службы; 

— разработка четырёхполюсной модели шестерённого насоса как эквивалентного источника колебаний, 
учитывающей приводную механическую и выходную гидравлическую системы; 

— разработка модели шестерённого насоса как эквивалентного источника колебаний, учитывающей его 
конструктивные особенности, процесс вытеснения, обратный гидроудар в межзубовую камеру, работу 
запертого объёма и перетечки из области нагнетания в область всасывания. 

При этом для подтверждения работоспособности предлагаемого подхода на весь класс машин объёмного 
вытеснения необходимо расширение номенклатуры испытуемых объёмных насосов. 
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